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Abstract 



We uncover an interesting phenomenon that neutrino flavor transformation in 
slowly varying matter density imitates almost exactly that of vacuum neutrino 
oscillation under suitably chosen experimental parameters. It allows us to 
have relatively large CP violating measure AP = ^(t'^ — > i^e) — P{i^fi ~^ 
Vg) which is essentially free from matter effect contamination. We utilize 
this phenomenon to design a low-energy long-baseline neutrino oscillation 
experiment to measure the leptonic CP violating phase. 
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Exciting discovery of neutrino oscillation in atmospheric neutrino observation [|I| and the 
persistent discrepancy between the observed and the calculated flux of solar neutrinos p| 
provide the strongest evidence for neutrino masses and lepton flavor mixing. Determination 
of all the mixing parameters, in particular the CP violating Kobayashi-Maskawa phase 0, 
is one of the most challenging goals in particle physics. 

This is the third in a series of works [^,^ in which we intend to explore possible (and 
hopefully experimentally feasible) ways of measuring leptonic CP violation in neutrino oscil- 
lation experiments. For early references and recent works on CP violation (or equivalently 
T violation), see for example, Refs. and Refs. HJ^, respectively. We focus in this paper 
the neutrino mass (difference) hierarchy suggested by the atmospheric neutrino observation 
and the MSW solutions [JlO[ of the solar neutrino problem ^j. In the standard three-flavor 
mixing scheme of neutrinos they exhaust all the independent mass difference squared: from 
atmospheric neutrino data, Am^g ^ Am^g = Arri^^^ ~ (2 — 5) x 10~'^ eV^ and from solar 
neutrino data, Amj^ = Am^^^^^ ~ (4-10) x IQ-*^ eV^ (SMA), (2 - 20) x IQ-^ eV^ (LMA) 
or (6 - 20) X 10"^ eV^ (LOW) where SMA, LMA, and LOW denote the small mixing an- 
gle, large mixing angle and the low Am^ MSW solutions, respectively |jll,12|. We use the 



notation in this paper as Am^j = mj — mj. 

Many authors including us examined the question of how to separate the genuine CP 
violating effect due to the leptonic Kobayashi-Maskawa phase from the fake one induced 
by matter effect In this paper we take a simple alternative strategy to look for a 

region of parameters in which the matter effect is "ignorable" in the first approximation. 
More precisely speaking, we will look for the solution of the question; is there region of 
tunable parameters in experiments, such as energy of neutrino beam, baseline length, etc. 
in which the neutrino oscillation probability, including its CP-odd term, are dominated by 
vacuum mixing effects? (See below on what we mean precisely by "dominated by vacuum 
mixing effects".) We will answer the question in the positive and finally end up with a 
proposal of experiment which utilizes low-energy neutrino beam of ~ 100 MeV and a 
megaton water Cherenkov detector to measure leptonic CP violation. 

We define the flavor mixing matrix U a.s = Uaii^i, where 1^0(0 = e,fi,T) and = 
1,2,3) stand for the gauge and the mass eigenstates, respectively. We take for convenience 
the representation of U as 



U 



Tse 



(1) 



C12C13 



C12C23 — Sl2'523'5l3e*'' -5230136*' 



■-512^23 "~ Ci2S23Si3e'" C12C23 — Si2S23Si3e'" 32301^6"" ^ (2) 

S12S23 - Ci2C23Si3e'" -C12S23 - Sl2C23Sl3e'" 0230136''' 

where Aj are SU(3) Gell-Mann's matrix and = diag{l, l,e*''). 

Let us first assume that the matter effect plays a minor role and consider the CP violating 
effect in vacuum. We will justify this assumption later. Under the mass difference hierarchy 
Ami3 ^ Am^2 which is implied by the solar and atmospheric neutrino data, the neutrino 
oscillation probability in vacuum can be written as 

'Al3L^^ tt* tt i ■ 2 f ^uL" 



2 • 2 

/33\ sm 



-2Jsin(Ai2L)[l - cos(Ai3L)] -h 4Jsin(Ai3L) sin^iAuL), 



(3) 
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where Ajj _ 
through J defined by 



One of the most significant feature of is that CP violation comes in 



J, 



af3;i,j 



lm[U^,U:^U;,U^,] (4) 
as it is the unique (in three-fiavor mixing scheme) measure for CP violation as first observed 



by Jarlskog |[T^ in the case of quark mixing. It takes the form in the parametrization we 



introduced above as J = ±ci2Si2C23S23C|3Si3 sin5, where the sign is positive for (e,/i) and 
(1, 2) and +(— ) corresponds to their (anti-) cyclic permutations of and 

We first observe that from the expression of J, if any one of the mixing angles is extremely 
small or very close to 7r/2 there is little hope in detecting the leptonic CP violation. For this 
reason, we will not deal with the case of SMA MSW solar neutrino solution. See, however, 
a remark at the end of this paper. 

We also notice that in order to have appreciable effect of CP violation. 



A12L = 0.26 



L 



10^5 / VlOOO km/ VlOO MeV 



E 



-1 



(5) 



should not be too small. If Am^2 



^"^solar Kiuch smaller than 10 ^ eV^ we have to 



either use a neutrino energy much lower than ~ 100 MeV which would be impractical in 
accelerator experiments, or take baseline length much longer than 1,000 km at the cost of 
lowering the beam fiux. For this reason, our discussion can accommodate neither the just-so 
vacuum oscillation solution nor the LOW MSW solution. Thus we assume in the following 
discussion that the nature chooses the LMA MSW solution to the solar neutrino problem. 

If we restrict ourselves into relatively short baseline, i.e., less than 1,000 km, one can 
safely assume that linear approximation sin (A12-L) — A12-L is valid. Therefore, in order that 
the CP violating (third) term in (^ is relatively larger than the 2nd term we must focus 
on the region where A12-L is smaller but not too smaller than unity. This consideration 
naturally leads us to the option of low-energy neutrino oscillation experiments with energy 
of order ~ 100 MeV. 

Next we must justify our assumption that the matter effect plays a minor role even for 
the CP violating effect. What we will show below is actually that whereas the Jarlskog factor 
and the energy eigenvalues are strongly modified by the matter effect it almost cancels out 
and does not show up in the observable quantities, the oscillation probabilities. 

We now discuss full system of three-flavor neutrino propagation in earth matter to under- 
stand the phenomenon of matter neutrino oscillation imitating vacuum oscillation. Toward 
this goal we develop an analytic framework based on perturbation theory under the adia- 
batic approximation We start from the neutrino evolution equation in the flavor 

basis: 
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(6) 



where a{x) = \/2G fN(,{x) indicates the index of refraction with Gp and A'^e being the Fermi 
constant and the electron number density, respectively. 
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We first note that at neutrino energy of ~ 100 MeV tliere is a liierarcliy among tlie 
relevant energy scales; 



2E 



The latter relationship holds because 
a(x) „ _ / p 



Am^/2E 



2.1 



> a 



2.7 g/cmV V0.5/ \10-^ eV 



2E 



Am? 



12 



E 



100 MeV 



Thanks to the mass hierarchy we can formulate the perturbation theory. 
We rewrite the evolution equation (H) as 



where u is defined by 



and Hq = diag{0, 0, A13) and 



J- A C 



a/3 



i^/3, 
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■=■12 


Cl2'5l2 


0" 




^13 





Cl3'5l3 


i/' = A12 


C12S12 


^12 





+ a{x) 
























_ C13S13 





'='13 



(7) 



(8) 



(9) 



(10) 



(11) 



The diagonalization of the 2x2 submatrix gives rise to the energy eigenvalues hi and the 
matter enhanced 612 as 



cl^aix) + A12 T V (cos 2012A12 - ac?3)2 + Af2 sin^ 2^12 



(12) 



sin2^f^ 



sin 26' 



12 



cos 2^12 -^c?3)2 + sin2 2^12 



(13) 



The matter "suppressed" for antineutrinos and the energy eigenvalues can be defined 
analogously by flipping the sign of a in (|12D and (piSl) . The resonance condition is satisfied 
in neutrino channel Q at 



Ci^a{x) = cos 26^12 A12 
which leads to the resonance energy of order of 100 MeV, 



(14) 



^ We note that one can take Aml2 always positive in the MSW mechanism as far as 612 is taken 
in its full range [0, 7r/2] |l2|. 
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Er = 48.6 COS 20 



P 



12 



V2.7g/cm' 



(15) 



for small value of the remaining vacuum mixing angle, ^^13 ^ 1. For the best fit values of the 
parameter of the LMA MSW solution, Amfg ~ (3.0 - 4.0) x 10^^ eV^ and sin^ 26^12 ~ 0.8 
]n[, we obtain ~ 60 - 80 MeV. 

The key to the matter enhanced 612 mechanism is the degeneracy of the zeroth order 
"energy" eigenvalue in the first 2x2 subspace of Hq. The degenerate perturbation theory 
dictates that one has to first diagonalize Hq + H' in this subspace to obtain the first-order 
corrected energy eigenvalues and the zeroth- order wave function. This last point is crucial. 
It gives rise to the matter enhanced CP violating effect that is free from the suppression by 
energy denominator. To our knowledge, this is the unique case of having CP violating effect 
which is not suppressed by any hierarchical ratios such as t^rn^^l 
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While our formalism is exactly the same as the one developed in Ref. |]14| the interpre- 
tation of the physical phenomena that occur is quite different from theirs. Also the nature 
of the resonance is quite different in large mixing angles. The resonance width can be esti- 
mated as ^ = tan 2^12 which means that ~ lEr ^ 120 MeV. Therefore, the resonances 
are so broad that they lose the identity as sharp resonances. 

We now show that despite the fact that 6'i2 could be strongly modified in matter, espe- 
cially for small vacuum ^12, the system mimics the vacuum neutrino oscillation even at the 
resonance where the effect of matter could be maximal. To understand this point we cal- 
culate the neutrino and antineutrino conversion probabilities P[y^l, and P[p^^ — > t'e)- 
Since the matter enhanced mixing angle 6'f| just replaces dyi in zeroth-order wave function it 
is straightforward to compute neutrino oscillation probabilities under the adiabatic approx- 
imation. They are nothing but the oscillation probabilities in vacuum but with 6'i2 replaced 
by 0*1, and Ajj by integrals over the energy eigenvalues, h\^2 and /is ~ A13. For example, 
the appearance probability P(z/^ — > z/g) reads 



4s23C?3S?3 sin2(^Ai3L) 



+C;^3 sin IQy^ (c23 — S23S;^3) sin 1^x2 + 2C23S23S13 COS b cos 26* 



12 



X sm 



-J (cos 26*12 



-2 JM(0i^, 5) sin J (cos 20 




A 



?.)2 + sin2 20i2Ai2l> 



-13 



12 



12 



A 



-c?3)2 + sin2 20i2Ai2L 
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(16) 



where Jm is the matter enhanced Jarlskog factor, JniOf^, S) = cos 0^ sin 0if C23S23C13S13 sin S, 
and we have averaged the rapidly oscillating piece driven by A13 in the CP violating term. 
The antineutrino transition probability -P(^'^ 



Ue) is given by the same expressions as 
above but replacing 0(2 by 0{2 and 6 by —6. 

Notice that at relatively short baseline, L < 1,000 km or so, the approximation sinx ~ x 
is valid. Then, the expressions of the oscillation probabilities approximately reduce to those 
in the vacuum because 



sin 20 




cos 2012 =F — cj.y + sin^ 20i2 A12 = sin 20i2 A 



A 



12- 



12 



(17) 



5 



Only mild a-dependence would remain due to the cos 5 term in (plq). Hence, as long as ^i2L 
is small, the oscillation probabilities of neutrinos and antineutrinos in matter imitate those 
in vacuum, independent of the mixing parameters and neutrino energy. 

We verify by numerical computations without using analytic expression in eq. (|16]) that 
the matter effect in fact cancels out. We tentatively take the values of the parameters as 
Am?3 = 3x10-3 eV^, sin^ 2^23 = 1-0, /^ml^ = 2.7x10"^ eV^, sin^ 2^12 = 0.79, sin^ 2^13 = 0.1 
and the Kobayashi-Maskawa phase 6 = tt/2. Our choice of 613 is within the newer CHOOZ 
bound Ijl^, but it may still be an optimistic one. With these parameters the Jarlskog factor 
in vacuum is given by J = 0.035, a small but non-negligible value. If we take 6 = — 7r/2 it 
corresponds, in a good approximation, to interchanging u and iy because the matter effect 
is only minor. 

In Fig. |1| we plot the oscillation probabilities for neutrino and anti-neutrino, and their 
difference, AP = Pliy^ Ve) ~ P{^pL ~^ ^e)i as a function of distance from the source with 
the neutrino energy E = MeV, for both z/^ — > v^. and z7^ — *• v^- It corresponds to the 
resonance energy in the neutrino channel. Insensitivity of the transition probability to the 
matter effect and how well it mimics the vacuum oscillation probability is clearly displayed 
in this figure. Although we sit on at the resonance energy of neutrino channel the features 
in the resonant neutrino flavor conversion cannot be traced in Fig. but rather we observe 
the one very much similar to the vacuum oscillation. We have also checked that the matter 
effect is small even if cos 5 term in ([T6|) is non-zero. The feature that it is a product of two 
harmonics with quite different frequencies can be understood by (^. We note that a large 
amplitude of | AP| (as large as ~ 0.3) in Fig. ^ can be qualitatively understood from the 
third term in our approximated analytic expression in eq. which implies that |AP| can 
be as large as 8J(~ 0.28) at A13L = vr. 

Now we discuss the possible experiments which utilize the imitating vacuum mechanism 
to measure leptonic CP violation. Measurement of CP violation at a few % level at neutrino 
oscillation experiments at ^ 100 MeV leaves practically the unique channel v^j, Ve- 
Experimentally, what we can actually measure is the number of events not the oscillation 
probability itself. Here we consider the ratio of the expected number of events due to z/^ — >■ 
and z7^ — >■ Vf. reactions in order to quantify the effect of CP violation. 

In such low energy appearance experiment we must circumvent the following difficulties; 
(1) smaller cross sections (2) lower flux due to larger beam opening angle, AO ~ 1 (E/lOO 
MeV)~^ radian. Therefore, we are invited to the idea of the baseline as short as possible, 
because the luminosity decreases as as baseline length grows. But of course, a baseline as 
long as possible is preferable to make the CP violating effect [the last term in (p!6|)] maximal. 
Thus we have to compromise. It appears to us that for E ~ 100 MeV, the baseline of ~ 
30-50 km seems preferable, as we will see below. 

Detection of low-energy neutrinos at a few % level accuracy requires supermassive de- 
tectors. Probably the best thinkable detection apparatus is the water Cherenkov detector 
of Superkamiokande type. Let us estimate the expected number of events at a megaton 
detector placed at L = 100 km. We assume that the neutrino beam flux 10 times as intense 



as (despite the difference in energy) that of the design luminosity in K2K experiment [19 



In the future it is expected that an 100 times more intense proton flux than KEK-ps seems 
possible at Japan Hadron Facility [^] but we assume only 10 times larger flux considering 
the fact that lower energy beam is more difficult to prepare. 
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The dominant i/p-induced reaction in water at around E = 100 MeV is not the famihar 



i6n ,,16/ 



z/g — e elastic scattering but the reaction on O, vl O — > e F ||18[. The cross section of the 

former reaction is about crfz/eC Uec) = 0.93 x 10^^^ ( ) cm^, while the latter is 

^ ^ VlOO MeV/ 

a(z/f O e-F) ~ lO^^^cm^ at E = 100 MeV |18|, which is a factor of 1000 times larger.[| 



But the number of oxygen in water is 1/10 of the number of electrons. So in net the number 

of events due to the reaction i^l^O e~F is larger than that of z/gC elastic scattering by 

a factor of 100. The number of ^^O in detector is given by 3.34 x 10^^ per kton of water. 

The neutrino flux at the detector located at L = 250 km is, by our assumption, 10 times 

more intense than the neutrino flux at Superkamiokande in K2K experiment. The latter is, 

/POT\ 

roughly speaking, 3 x 10^ \~[^o') where POT stands for proton on target. Therefore, 

the expected number of events N assuming 100 % conversion of to z/g is given by 



N ~ 6300 



)""(nLi)(w)- 



100 km 



In the antineutrino channel, the dominant reaction is UeP —>■ e~^n with cross section 
cr ~ 0.4 X 10~^^cm^ at E = 100 MeV. The event number due to this reaction, assuming the 
same flux of t'^ as z/^, is similar to that of (|18]) because the cross section is about half but 
there are two free protons per one oxygen. There is additional oxygen reaction I'l^O — >■ e'^N 
with approximately factor 3 smaller cross section than that of i^l^O |18| . 



In order to estimate the optimal distance, we compute the expected number of events 
in neutrino and anti-neutrino channels as well as their ratios as a function of distance with 
taking into account of neutrino beam energy spread. For definiteness, we assume that the 
average energy of neutrino beam {E) = 100 MeV and beam energy spread of Gaussian type 
with width aE = 10 MeV. We present our results in Fig. 0. 

If there is no matter and CP violating phase effects the ratio R = Nlu^ z/e)/A^(z/^ — »• 
z/g) must not vary with the distance but must be a constant which is (j)u{<^u) / <t>p{'^p) where 
(pu and ai, (cxp) are total flux and cross section averaged over neutrino (anti-neutrino) 
energy, respectively. Taking the error (only statistical one for simplicity) into account, L 
must not be too short or too long. The optimal distance turned out to be rather short, 
L ~ 30 — 50 km, much shorter than that of the current long-baseline experiments, (see Fig. 
2.) At such distance the significance of the CP violating signal is as large as 3-4 cr. 

Roughly speaking, the deviation of R from the case without CP violating effect, AR, 
is proportional to AP but what is relevant is the quantity AP x (A^ is the expected 
number of events in the absence of CP violating effect) which signals if the CP violating 
effect is significant or not. From eq. (^ we can see that when the distance is small so that 
A13L <^ 1, the terms proportinal to J (or AP) behaves like AP oc [L/E)^. This implies 
that assuming (p^, oc l/L"^, AP x N cc APa{E)/L'^ oc a{E)L/E^ is an increasing function of 



^ Notice that the energy dependence of I'l^O reaction is very steep below E = 100 MeV. For 
example, it is smaller by an order of magnitude at S = 60 MeV. Because the energy dependence 
of AP is mild we are invited to take larger energy than the resonance energy. 
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L. When the distance is larger and A13L is not small, AP oc (L/E) and this implies that 
AP X N (X a{E)E/L is a decreasing function of L. 

Hence, the significance of CP violating effect is expected to be maximal at some distance 
L where AP x takes the largest value. In our case it corresponds to L ~ 30 — 50 km for 
the mixing parameters and neutrino energy we consider, as indicated in Fig. 2. 

Let us also make some comments on the validity of our choice of low neutrino energy (as 
low as ~ 100 MeV) from the view point of the expected number of events. Suppose that 
the relevant cross section (t{E) behaves like (t{E) oc -E" and we can neglect matter effect. 
In this case, is ~ f{L/E) X where /(L/E) is some function of L/E. From this we 

can conclude that if a = 2 the same number of events can be expected as long as the value 
of L/E is the same if the same intensity of neutrino beam is assumed because the increase 
in the cross section with higher energy can be compensated by the decrese of the flux due 
to larger distance. On the other hand, if a 7^ 2 such scaling law does not apply. If a > 2 
{a < 2) the one can expect larger number of events with higher (lower) energy if the value 
of L/E is kept constant. For neutrino energy above 100 MeV but below 1 GeV, for neutrno 
channel, roughly speaking a ~ 2 whereas for anti neutrino channel a ~ 1 P2[. Therefore, 



our choice of lower energy seems to be preferable. For energy well above GeV, a ~ 1 for 
both channels and lower energy is preferable. 

We conclude that the experiment is quite feasible under such intense neutrino beam and 
a megaton detector. Fortunately, the possibility of constructing a megaton water Cherenkov 
detector is already discussed by the experimentalists pl |. 

Finally, a few remarks are in order: 

(1) Imitating vacuum mechanism does work also for the SMA solar MSW solution after 
averaging over the rapid oscillations due to larger Am^. Because of the tiny vacuum angle, 
however, the number of events is smaller by a factor of ~ 2500 than the LMA case and the 
experiment does not appear feasible with the same apparatus. 

(2) An alternative way of measuring CP violation is the multiple detector method which 
may be inevitable if either one of or z/ beam is difficult to prepare. It utilizes the fact that 
the first, the second, and the third terms in the oscillation probability ([T6|) have different L 
dependences, ~ L- independent (after averaging over energy spread of the neutrino beam), 
~ L^, and ~ L, respectively, in the linear approximation.^ 

(3) Intense neutrino beams from muon storage ring at low energies, proposed as PRISM, 



Phase- Rotation Intense Secondary Mesons |23], would be an ideal source for neutrinos for 



the experiment proposed in this paper. Of course, it requires identification of z7g from i/g by 
some methods, e.g., by adding Chlorine (^^Cl) into the Water Cherenkov detector to make 
it sensitive to the characteristic ~ 8 MeV 7 rays arising from the absorption of neutron into 
the Chlorine followed by UeP — ^ e^n reaction. 



^ The minimal two detector methods Q is most powerful if the far (near) detector is placed 
at A12 = ^(f ) because then the difference of the CP violating term in ( |l^ ) between the far 
and the near detectors is maximal. However, it requires the basehne length of far detector as 



L = 6, 300 n km which may or may not be too large to realize. 

\3xlO-5eV2/ V 100 MeV/ ^ ^ ^ 
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(4) We would like to urge experimentalists to think more about the better supermassive 
detection apparatus than water Cherenkov for highly efficient and accurate measurement of 
low energy neutrinos. 
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FIG. 1. Oscillation probability for (a) neutrinos, -P(f^ I'e)^ (b) anti-neutrinos, P{i>fji i^e), 
and (c) their difference, AP(f^ v^) = P{Vfj, i'e) — P{v'n u^.) with fixed neutrino energy Ey = 
60 MeV, are plotted as a function of distance from the source. The mixing parameters are fixed to 
be Amfa = 3 x IQ-^ eV^, sin^ 2^23 = 1-0, Amfg = 2.7 x IQ-^ eV^, sin^ 2^12 = 0.79, sin^ 2^13 = 0.1 
and 6 = 7r/2. We take the matter density as p = 2.72 g/cm^ and the electron fraction as = 0.5. 
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FIG. 2. Expected number of events for (a) neutrinos, Nlu^ 
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(b) anti-neutrinos, 



Ve), and (c) their ratio R = N[u^ — > Ue)/N[u^ u^) with a Gaussian type neutrino 



energy beam with {E^) = 100 MeV with cr = 10 MeV are plotted as a function of distance from 
the source. Neutrion fluxes are assumed to vary as ~ 1/L^ in all the distance range we consider. 
The mixing parameters as well as the electron number density are fixed to be the same as in Fig. 
1. The error bars are only statistical. 
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